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Abstract 

Characteristics of hadron production in diffractive deep-inelastic positron-proton 
scattering are studied using data collected in 1994 by the HI experiment at HERA. 
The following distributions are measured in the centre-of-mass frame of the pho- 
ton dissociation system: the hadronic energy flow, the Feynman-x (xp) variable 
for charged particles, the squared transverse momentum of charged particles (p?p), 
and the mean as a function of xp. These distributions are compared with re- 
sults in the j*p centre-of-mass frame from inclusive deep-inelastic scattering in the 
fixed-target experiment EMC, and also with the predictions of several Monte Carlo 
calculations. The data are consistent with a picture in which the partonic structure 
of the diffractive exchange is dominated at low Q 2 by hard gluons. 
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1 Introduction 



Studies of the 1992 deep-inelastic electron-proton scattering events (DIS) at HERA |jj revealed 
the presence of "rapidity-gap events" - events of the form ep — > eXY in which the hadronic 
final state consists of two parts, X and Y, separated by a large region in pseudorapidity in 
which no hadrons are observed. This is illustrated in figure [I]. The masses Mx and My of these 
two systems, separated by the largest rapidity gap in the event, are thus small compared to 
W, the invariant mass of the "f*p system. The system Y in these events consists of a proton 
or other low-mass hadronic state and has a momentum similar to that of the incoming proton. 
The magnitude of the square of the 4-momentum transferred to X from the proton, is 
small lGeV 2 ). The contribution of these events to the ep interaction cross section has 
been measured in terms of a diffractive structure function, F^ 3 ', using the HERA data of 
1993 [||. The more precise measurements made with the 1994 data |3| indicate that 

the cross section for rapidity-gap events may be parameterised as a diffractive contribution, 
from pomeron (JP) exchange, together with a contribution from meson exchange. The pomeron 
contribution dominates for small values (< 0.05) of the variable xjp, which is the fraction of the 
incoming proton's longitudinal momentum carried by the exchange. Furthermore, the pomeron 
may be interpreted as having partonic structure. 

A QCD study of parton distribution functions ||], evolved according to the DGLAP equa- 
tions reveals the preference of the F®^ 1 data for a pomeron that is dominated by a "hard- 
gluon" parton distribution at the starting scale of Qq = 3GeV 2 (fits 2 and 3 in ||), i.e. a 
distribution with a large contribution from gluons carrying a significant fraction of the momen- 
tum of the pomeron. A pomeron model with only quarks at Qq (fit 1 in Q) fails. 




Figure 1: Schematic picture of a deep-inelastic interaction ep — > eXY with a "rapidity gap" devoid of hadronic 
energy. The high-Q 2 virtual photon interacts with a colourless space-like entity of squared 4-momentum t, whose 
longitudinal momentum as a fraction of that of the target proton is xp; the direction of the system Y is close to 
that of the target proton. 

The studies of the diffractive final state presented here may be regarded as specific tests of the 
pomeron structure extracted from F% , complementary to other final-state analyses 0, 0, M. 
A gluon-dominated pomeron is expected to interact largely by boson-gluon fusion (BGF); the 
transverse momenta of the outgoing partons relative to the photon direction in the 7*iP centre- 
of-mass (CM) frame (i.e. the CM frame of the photon dissociation system X), as well as the 
amount of gluon radiation and consequently the parton multiplicity, are therefore likely to be 
greater than in the case of a quark-dominated pomeron. 



The partonic structure of the pomeron should thus be reflected in the energy-flow and 
charged-particle distributions in the 7*iP CM frame. These distributions are compared here 
with data on deep- inelastic \xp interactions in the 7*p CM frame at the EMC experiment. This 
provides a comparison of the structure of the pomeron with that of the proton. 

The distributions are also compared with predictions from a calculation using a pomeron 
with parton distributions from the best QCD fit to F% (fit 3). Comparison with predictions 
from a calculation with only quarks in the pomeron at the starting scale of the QCD fit (fit 1) 
demonstrates the sensitivity of the measurements to the parton distributions of the pomeron. 

Another approach to diffractive DIS is provided by the photon dissociation picture ||. In 
the rest frame of the proton, the photon fluctuates a long time before the interaction into a Fock 
state with definite parton content (qq, qqg, . . . ), and this partonic state scatters diffractively off 
the proton. Scattering from a quark in the pomeron corresponds to an Aligned Jet Model H 
topology in this picture, whereas scattering from a gluon corresponds to Fock states with one 
or more gluons and results in less pronounced alignment. 

A non-diffractive model for the production of events with a large rapidity gap has also been 
proposed, using a conventional picture of DIS based on scattering from a single parton within the 
proton, followed by soft colour interactions [1C, 11|. The data are also compared with predictions 
of this model. 



2 Detector, Event Selection and Kinematics 

The data used here were collected in 1994 with the HI detector at the HERA collider, which 
operated with 27.5 GeV positron and 820 GeV proton beams. The HI detector is described in 
detail elsewhere [jl^]; those components of importance for the analyses presented here are briefly 
mentioned in the following. The laboratory coordinate system has its origin at the nominal 
interaction point and its z axis in the proton beam direction, also called the forward direction. 
The pseudorapidity is defined as r\ = — In tan |, where 9 is the angle with respect to the z axis. 

In the forward and central regions the interaction point is surrounded by a system of tracking 
detectors - interleaved drift and multi-wire proportional chambers - which cover the pseudo- 
rapidity range —1.5 < rj < 2.8 and the full azimuth. The momenta of charged particles are 
determined for this analysis from their track curvature in the central jet chamber (CJC) in the 
uniform magnetic field of strength 1.15 T generated by a superconducting solenoid. The CJC has 
a resolution of cr pT /px ~ O.OOQpx © 0.008 (where px is in GeV and the constant term describes 
the contribution from multiple scattering at pt = 0.5 GeV) and ag ~ 20 mrad. The solenoid sur- 
rounds the trackers and the fine-grained liquid-argon (LAr) calorimeter, which covers the range 
— 1.5 < i] < 3.6 and is used to measure energies in the hadronic final state with a resolution of 
<Je/E ~ 0.5/\/^©0.02 (where E is in GeV) [O]. Charged-particle detection in the backward 
region, —3.1 < 77 < —1.5, is provided by the backward proportional chamber (BPC). Behind 
this, the backward electromagnetic calorimeter (BEMC) completes the calorimetric coverage in 
the range —3.5 < r) < —1.5. 

Several subdetectors in the forward region are used in this analysis to tag particles emitted 
close to the proton direction. These detectors cover a range larger than their purely geometrical 
acceptance, due to the effect of secondaries resulting from the scattering of primary particles 
in the beam pipe and adjacent material. The copper-silicon plug calorimeter covers the range 
3.5 < rj < 5.5. The first three double layers of drift chambers in the forward muon detector 
(FMD) are sensitive to particles in the range 5.0 < r\ < 6.5. The proton-remnant tagger (PRT), 
comprising a set of double layers of scintillators situated around the proton beam pipe, covers 
the range 6.0 < rj < 7.5. 



The luminosity is determined from the rate of Bethe-Heitler interactions detected in a mon- 
itor downstream of the main detector in the positron direction. The total luminosity analysed 
here amounts to approximately 2pb _1 . 

All events included in this analysis are required to have a scattered positron reconstructed 
with an energy greater than 12GeV, emitted in the angular range 156° < 9 < 173°. The 
scattered positron is identified by looking for the highest-energy cluster in the BEMC. This 
cluster is required to be associated with a hit in the BPC and to pass a cut on its transverse size. 
The trigger used to select these events for read-out requires the presence of an energy deposit of 
at least 4 GeV in the BEMC. The efficiency of the BEMC trigger for positron energies larger than 
10 GeV is known to exceed 99% |13|| . For events passing the above selection the photoproduction 
background is less than 1%. 

The data are selected to lie in the kinematic region 7.5 < Q 2 < 100 GeV 2 and 0.05 < y < 0.6. 
These variables are defined as Q 2 = —q 2 and y = (q ■ p)/{k ■ p) ~ (W 2 + Q 2 )/s, where p, k and 
q are the 4-momenta of the incoming proton, positron and photon respectively, W is the CM 
energy of the 7*p system, and s is the squared CM energy of the ep system. The variables Q 2 
and y are calculated from the measured energy and angle of the scattered positron. 

Diffractive events are recognised by the rapidity gap between the outgoing hadronic systems 
X and Y . This gap is tagged by the absence of signals above noise in the forward components 
of the HI detector |3|, [14|]: the plug calorimeter, the FMD and the PRT. There must also be no 
energy deposit of more than 400 MeV in the forward region [r] > 3.0) of the LAr calorimeter. 

The variable xp is defined as 

q-P 

Xp = 

q-p 

where P is the 4-momentum exchanged between systems X and Y, and is reconstructed using 
the relation 

x P ~ — (E + Pz)- 
itj v e+X 

Here, E p is the energy of the incoming proton, E and p z are the energy and longitudinal 
momentum of each final-state particle in the laboratory frame, and the sum runs over the 
scattered positron e and all detected particles in the photon dissociation system X. The particles 
are reconstructed using a combination of tracks and calorimeter clusters, with an algorithm for 



track-cluster association avoiding double counting 15]. In order to enhance the contribution 
from pomeron exchange and to reduce the contribution from meson exchange, the requirement 
xp < 0.025 is made. The fraction of the exchanged momentum P carried by the struck quark 
is given by (3, where 

Q 2 Q 2 



^ 2q-P Q 2 + M\' 

Neither the squared momentum transfer t nor the mass My of the hadronic system Y is 
measured here. However, the requirement of the absence of activity in the forward detectors 
imposes the approximate restrictions My < 1.6 GeV and \t\ < lGeV 2 , and the results are 
corrected to this kinematic region. Since My is not measured, it is not possible to distinguish 
events containing an elastically scattered proton from those in which the proton dissociates into 
a low-mass state. 

The analyses are performed in the 7*iP CM frame, the 4-momentum of the 7*iP system 
being reconstructed as V = q + xpp [16|. The transverse momentum of the pomeron with 



respect to the direction of the incoming proton is not measured, but can be neglected since this 
has no significant effect on the measurements presented here. 



The energy-flow distribution is measured using energy deposits reconstructed from clusters 
of cells in the LAr and BEMC. The clusters are required to have an energy in the range 
< E < 25 GeV and to lie in the angular range —3 < rj < 4. The charged-particle distributions 
are measured using only particles that are detected in the CJC, originate from the ep interac- 
tion vertex in the angular region —1.31 < 77 < 1.31, and have a transverse momentum in the 
laboratory frame of 0.15 < px < 10 GeV, but the results are fully corrected with Monte Carlo 
simulations (see below) to cover the whole hadronic system. 

Two variables are calculated for each charged particle, both defined in the CM frame of the 
7*iP system. The Feynman-x variable, xjr, is defined as 

2pj! 

where pjj is the component of the particle's momentum in the direction of motion of the incoming 
photon. In non-diffractive DIS, xp is defined in the j*p CM frame, and the denominator is W. 
Positive xf corresponds to the current hemisphere and negative xf to the target fragmentation 
hemisphere. The other variable used here, p T , is the transverse momentum of the particle with 
respect to the photon direction. 



3 Monte Carlo Models for Diffractive Interactions 



Monte Carlo calculations employing the event generator RAPGAP 2.02 [17] are used in conjunc- 
tion with a detailed simulation of the HI apparatus to correct the experimental distributions for 
the acceptance and resolution of the detector, and later to compare the results with theoretical 
expectations. The data are also compared with the predictions of the Monte Carlo generator 



LEPTO 6.5 |Tl| 



The RAPGAP model treats diffractive interactions as inelastic eiP collisions, the pomeron 
being modelled as an object with partonic substructure. Scattering on mesons is also included. 

For scattering from quarks, the lowest-order diagram considered by RAPGAP at the parton 
level is 0(a) eq — > eq scattering (figure Qa); a higher-order, 0(aa s ), process is QCD Compton 
(QCD-C) scattering eq — > eqg (figure |2]b). For scattering from gluons, the lowest-order 0(aa s ) 
process is boson-gluon fusion eg — > eqq (figure ||c) . A cut p\ > 2 GeV 2 is applied in order to 
avoid divergences in the 0(aa s ) matrix elements for massless quarks, px being the transverse 
momentum of the outgoing partons with respect to the photon direction in the CM frame of the 
hard subprocess. For p T < 2 GeV 2 , only eq — > eq scattering is used. In the RAPGAP model, a 
"pomeron remnant" is implicit at the parton level, consisting of a quark (antiquark) in eq (eq) 
scattering and a gluon in the eg case (see figure . The renormalisation and factorisation scale 
fj, 2 is set to Q 2 . Higher-order corrections at the parton level are treated using leading- log parton 
showers (PS) JTH] as implemented in RAPGAP. Subsequent hadronisation is simulated according 



to the Lund string model in JETSET [19|. The dependence of the acceptance corrections on 
the method used to handle higher-order corrections in RAPGAP is studied with a separate 
calculation in which QCD-C and higher-order processes are simulated by the colour-dipole (CD) 
model |2(J. RAPGAP is interfaced to HERACLES for the simulation of QED radiative 
effects. 

For the comparisons presented here, two sets of parton distributions for the pomeron are used 
in the RAPGAP generator. The first is taken from the best fit to the Q 2 and (5 dependence of the 
diffractive structure function (fit 3 in ||). It is characterised by a hard gluon distribution 

at the starting scale for DGLAP evolution, Qq = 3 GeV 2 , in which gluons carry > 80% of the 




Figure 2: Elementary processes included in the RAPGAP simulation of ep diffractive interactions: (a) lowest- 
order, 0(a), e scattering from a quark, leading to a two-parton configuration; (b) one of the diagrams for the 
QCD-C process, 0(aa„); (c) lowest-order, 0(aa s ), e scattering from a gluon by BGF, leading to a three-parton 
configuration. 



total momentum of the diffractive exchange^} The second is taken from a fit in which only 
quarks are permitted to contribute to the partonic structure of the pomeron at Q 2 = Qq (fit 1 
in §). The latter choice does not provide a satisfactory description of the measurements; 
it is used here to demonstrate the sensitivity of the results to the parton distributions of the 
pomeron. The predictions of RAPGAP with these two sets of parton distributions are hereafter 
referred to as RG-F?? (fit 3) and RG-F?? (fit 1) respectively. The meson structure function is 
taken to be that of the pion p2[ . In both cases, the parton distribution functions evolve with 
Q 2 according to the DGLAP equations. The partons are treated as having no intrinsic k,T- 

In the LEPTO 6.5 model, rapidity-gap events are generated by the soft colour interaction 
(SCI) mechanism. Deep-inelastic scattering is modelled using the matrix elements for processes 
up to 0(aa s ), as in RAPGAP, but with partons coming directly from the proton according to 
the MRS(H) parameterisation f23[| of the proton structure function, in which the renormalisation 
and factorisation scale fj, 2 is again set to Q 2 . The divergences in the 0{aa s ) matrix elements are 
avoided using the cuts s > s m i n and z q < z q ^ m i n < 1 — z q , where s is the CM energy of the hard 
subprocess, z q = (p • p q )/(p ■ q) and p q is the 4-momentum of one of the outgoing partons. The 
parameters s m i n and z q ^ m i n are set to 4GeV 2 and 0.04 respectively. Higher-order corrections 
are treated with the PS method, and hadronisation follows the Lund string model. Further 
non-perturbative interactions take place as the outgoing partons pass through the colour field 
of the proton. These soft colour interactions can result in a hadronic final state comprising two 
colour-singlet systems separated by a rapidity gap. 



4 Acceptance Corrections and Systematic Errors 



The data are corrected for the acceptance and resolution of the HI apparatus using events 
generated by RAPGAP 2.02 17] with a hard-gluon pomeron structure function. The following 
sources of systematic error are taken into account; the errors shown in brackets are those on the 
energy flow and the charged particle d istributions respectively: 

1 Fits 2 and 3 in || give very similar predictions for the distributions studied here. Only the predictions from 
fit 3 are shown in this paper. 



— an uncertainty in the LAr calibration for hadronic energy of ±4%, leading to an average 
error in the experimental distributions of ±(6%, 2%); 

— an uncertainty in the BEMC calibration for hadronic energy of ±20%, leading to an average 
error of ±(6%, 0.7%); 

— an uncertainty of ±1% in the positron energy scale, leading to an average error of 
±(3%, 0.5%); 

— an uncertainty of ±lmrad in the positron polar angle, leading to an average error of 
±(0.5%, 0.5%); 

— an uncertainty in the xjp dependence in the Monte Carlo used to calculate the acceptance 
corrections, taken into account by reweighting the generated events with the function x^' 2 , 
resulting in an average error of ±(8%, 3%); 

— an uncertainty in the (3 dependence in the Monte Carlo, taken into account by reweighting 
the generated events with the function (a -1 — a) (3 + a, where a has the range 0.5 to 2, 
resulting in an average error of ±(6%, 1%); 

— an uncertainty in the t dependence in the Monte Carlo, taken into account by reweighting 
the generated events with the function e ±2 *, resulting in an average error of ±(5%, 2%); 

— uncertainties arising from the method used to treat higher-order corrections in the Monte 
Carlo; this results in an error of ±(10%, 10%), evaluated by taking the full difference 
between the results obtained using the PS and CD models as an estimation of the extra- 
polation uncertainty from the Monte Carlo. 



In the charged-particle analysis, the track selection criteria were varied to allow for the imperfect 
description of the Central Tracker in the Monte Carlo. The selection was varied by increasing 
the angular range of the tracks to —1.50 < 77 < 1.64 and by requiring the radial track length 
to be greater than 15 cm and the number of hits on a track N^ts to be greater than 10. This 
results in an average error of ±5%, rising in the low-xp region to 20% in the xp distribution 
and 15% in {p^) in the "seagull plot" (figure ||). 



5 Results 

Resultsn are shown for the kinematic range 7.5 < Q 2 < 100 GeV 2 , 0.05 < y < 0.6, xjp < 0.025, 
\t\ < lGeV 2 and My < 1.6 GeV. All distributions are fully corrected for the effects of the 
acceptance and resolution of the HI apparatus, and statistical and systematic errors on the data 
points are combined quadratically. The inner vertical error bars indicate the statistical error 
and the outer bars the total error. Due to the xjp cut, the meson-exchange contribution is small 
(< 7% in the RG-Fj? (fit 3) model) and does not affect the conclusions. 

It is shown below that, in Monte Carlo calculations, the characteristic features of the parton 
content of the pomeron are reflected in the distributions of various final-state observables. The 
data are therefore compared with data on ~f*p interactions at relatively high Bjorken-x values 
(x > 0.01), where the structure of the proton is dominated by quarks. The data used for this 
comparison are taken from \ip interactions at the EMC experiment [0] , with the mean j*p CM 
energy, (W) = 14 GeV, of the EMC inclusive DIS data similar to the mean 7*iP CM energy, 
(M x ) « 12 GeV, of the HI diffractive data. 

2 The data are available in numerical form on request and have also been submitted to the Durham HEPDATA 



database http : //durpdg . dur . ac . uk/HEPDATA 



The experimental data are also compared with the predictions of the RAPGAP model 
with the two different pomeron structures from the DGLAP fits to F®^ 1 - the preferred 
hard-gluon structure (RG-F^ (fit 3)) and the structure with only quarks at the starting scale 
(RG-F^ (fit 1)) - as well as with the soft colour interaction (SCI) model as implemented in 
LEPTO 6.5. 

In the CM frame of the 7*iP system, a model with a quark-dominated pomeron results 
in events in which the struck quark and the pomeron remnant tend to be close to the 7*iP 
axis. Large contributions to the transverse momentum can arise from QCD corrections such 
as QCD Compton scattering, but the rate of these is suppressed by a factor a s . In contrast, 
if the pomeron is a gluon-dominated object, more transverse momentum and energy flow are 
produced at lowest order. This arises because the quark propagator in the BGF process can 
have non-zero virtuality, so the quark-antiquark pair from the hard subprocess is not necessarily 
aligned along the 7*iP axis. Evidence for this effect is seen in a recent analysis, using the same 
data, of the distribution of the pr of thrust jets relative to the j*1P axis ||. In both the quark- 
and gluon-dominated cases, small contributions to the transverse momentum are generated by 
the intrinsic kx of the partons, by the hadronisation process, and by particle decays. However, 
it is known from the average thrust distribution shown in || that a quark-dominated exchange, 
even with large intrinsic kx, would not be able to explain the relatively low value of the average 
thrust. 

Figure ^ shows the event-normalised energy flow 1/N dE/dr/* in three different regions of 
Mx, rj* being the pseudorapidity relative to the direction of motion of the incoming photon 
in the CM frame of the j*1P system and iV being the number of events. The distribution is 
approximately symmetrical about rj* = 0, with similar levels of energy flow in the two hemi- 
spheres. At higher masses, Mx > 8 GeV, a two-peaked structure is seen, indicating that the 
major topological property of these events is of a 2-jet nature (current jet and remnant jet). 

The data are well described by the RG-F!? (fit 3) modef], and for M x > 8 GeV by 
LEPTO 6.5, whereas the RG-F^ (fit 1) model predicts too much energy flow at the largest 
accessible pseudorapidity and fails to account for the observed energy flow in the central region, 
rj* « 0. 

These features can be related to the gluon content of the pomeron and to the role of BGF. 
The transverse momentum generated by this process reduces the accessible pseudorapidity range 
along the j*1P axis at fixed Mx and induces additional energy flow in more central regions. In 
addition, more energy is expected in the central region for a gluon-dominated than for a quark- 
dominated object because of the exchange of a gluon in BGF; this results in a colour octet-octet 
field between the outgoing quarks and the gluonic remnant (which all together form an overall 
colour singlet), giving rise to enhanced soft gluon radiation compared to the triplet-antitriplet 
field between the struck quark and the remnant in eq scattering ]l6|, ^5|. This effect has been 
seen for gluon jets produced in e + e~ interactions |26|]. This enhancement is also apparent in the 
virtual photon dissociation picture where the process corresponding to BGF involves a Fock 
state containing a low-energy gluon. 

The production of charged particles is studied using the distributions of the variables xp and 
Pj?, defined in section [2|. The results are shown for the restricted mass range 8 < Mx < 18 GeV 
in order to make the data more directly comparable with the EMC j*p results. 

The xf distribution, normalised by the number of events N, is shown in figure |I[ A feature 
of the diffractive data is the similarity between the xf distributions in the two hemispheres, in 
contrast to the strong asymmetry of the *y*p data. The asymmetry in the latter case is explained 

3 The Pt cut in the generator was varied from 2 up to 4 GeV 2 without affecting the conclusions drawn from 
the distributions presented in this paper. 




Fi gure 3: The energy-flow distribution in the 7*JP CM frame for three different Mx intervals, in the kinematic 
region 7.5 < Q 2 < 100 GeV 2 , 0.05 < y < 0.6, x F < 0.025, \t\ < 1 GeV 2 and M Y < 1.6 GeV. Positive rj* 
corresponds to the direction of motion of the incoming photon. The statistical error bars are hidden by the 
symbols. The points are plotted at the centre of the bin in the horizontal coordinate, and the horizontal bars 
indicate the width of the bin. Also shown are the predictions of the following Monte Carlo models: RAPGAP with 
the hard-gluon pomeron structure taken from fit 3 in |] (RG-F? (fit 3)); RAPGAP with the pomeron structure 
containing only quarks at the starting scale, taken from fit 1 in J3j (RG-F? (fit 1)); and the soft colour interaction 
model as implemented in LEPTO 6.5. The particles are taken as massless for this distribution. 



by the requirement of baryon number conservation in ep interactions. This results in a large 
contribution from baryons in the region xf ~ — 0.4 J24| |, 

The xf distribution is well described by the RG-FJ? (fit 3) and LEPTO 6.5 models, whereas 
the RG-F?? (fit 1) calculation fails in that it predicts too little particle production in the central 
region [xf ~ 0) and too much particle production at large \xf\- These features are related to 
those observed in the energy-flow distribution: longitudinal phase space is restricted because of 
the transverse momentum generated by the BGF process, and there is additional radiation in 
the central region for a gluon-dominated object. 

The p^p distribution, normalised by the number of events N, is shown in figure ||. Results 
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Figure 4: The Feynman-a; (xf) distribution, showing HI 7*_P data in the -y*P CM frame in the kinematic 
region 7.5 < Q 2 < 100 GeV 2 , 0.05 < y < 0.6, x P < 0.025, 8 < M x < 18 GeV, \t\ < 1 GeV 2 and My < 1.6 GeV, 
together with EMC fip DIS data in the 7*p CM frame. Positive xf corresponds to the direction of motion of 
the incoming photon. Each point is plotted at the horizontal position where the predicted probability density 
function has its mean value in the bin, as prescribed in | p7| . The uncertainty in the horizontal position due 
to the model-dependence of the shape of the distribution is indicated by horizontal bars, but in most cases is 
contained within the width of the symbol. The bin boundaries are marked by vertical dotted lines. Also shown 
are the predictions of the following Monte Carlo models for the HI data: RAPGAP with the hard-gluon pomeron 
structure taken from fit 3 in || (RG-F!? (fit 3)); RAPGAP with the pomeron structure containing only quarks at 
the starting scale, taken from fit 1 in || (RG-F° (fit 1)); and the soft colour interaction model as implemented 
in LEPTO 6.5. 



are shown for the range 0.2 < xf < 0.4. This restricts the data to the current region, where a 
comparison between HI and EMC data is more meaningful, and matches the xp range chosen 
by EMC. The number of high-p^ charged particles (p^ > 2.0 GeV 2 ) is significantly higher in the 
7*iP data than in the EMC ^y*p data. This points to a larger contribution from scattering from 
gluons in the diffr active case than in inclusive DIS. 

The py? distribution is well described by the RG-F?? (fit 3) calculation, which predicts a 
harder spectrum than the RG-F?? (fit 1) model because of the role of the BGF process. The 
data are also well described by LEPTO 6.5. The RG-F?? (fit 1) model does not describe the 
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Figure 5: The p?j? distribution, showing HI 7*_P data in the -y*P CM frame in the kinematic region 7.5 < Q 2 < 
100 GeV 2 , 0.05 < y < 0.6, x P < 0.025, 8 < Mx < 18GeV, \t\ < lGeV 2 , M Y < 1.6 GeV, and 0.2 < x F < 0.4, 
together with EMC [ip DIS data in the 7*p CM frame. Each point is plotted at the horizontal position where 
the predicted probability density function has its mean value in the bin, as prescribed in | p7| . The uncertainty 
in the horizontal position due to the model-dependence of the shape of the distribution is indicated by horizontal 
bars, but in most cases is contained within the width of the symbol. The bin boundaries are marked by vertical 
dotted lines. Also shown are the predictions of the following Monte Carlo models for the HI data: RAPGAP with 
the hard-gluon pomeron structure taken from fit 3 in |] (RG-F? (fit 3)); RAPGAP with the pomeron structure 
containing only quarks at the starting scale, taken from fit 1 in J3j (RG-FJ? (fit 1)); and the soft colour interaction 
model as implemented in LEPTO 6.5. 



distribution well at high p^ 2 . It should also be noted that the RG-F^ (fit 1) predictions for the 
p?P distribution are in good agreement with the EMC 7*p data. This supports the argument 
that the EMC data (at a mean Q 2 of 12 GeV 2 ) are dominated by quarks and that the higher 
in the diffractive data is not an effect of the larger Q 2 at HI (where the mean Q 2 in this 
analysis is 25 GeV 2 ). 

Confirmation of the general trends discussed above is obtained from the "seagull plot" , shown 
in figure |(| in which the mean transverse momentum squared, (p^ 2 ), is plotted as a function of 
xp. It is observed that (pr) is significantly higher in the HI diffractive data than in the EMC 



j*p data at all but the highest xp values, confirming the conclusions drawn from the comparisons 
in figure [|. The seagull plot also shows a greater degree of symmetry than is present in the j*p 
data (cf. figure ||). As well as the effect of baryon number conservation, this is indicative of 
more parton radiation in the target fragmentation region. This is consistent with having a more 
point-like partonic system in the target fragmentation hemisphere in the diffractive data than in 
the j*p data, where the extended nature of the proton remnant leads to a restricted phase space 
for parton radiation. Note, however, that there might be an indication of a small asymmetry 
between the two hemispheres in the energy- flow distribution (figure ||). 

The significant increase in (pp 2 ) as \xp\ increases from to ~ 0.5 is well described by the 
RG-FJ? (fit 3) calculation and reasonably well described by LEPTO 6.5, whereas the (p?) values 
predicted by the RG-F?? (fit 1) model are too low across most of the xf range, in accordance 
with the above discussion of figure || 

6 Summary and Conclusions 

Energy-flow and charged-particle spectra have been measured in diffractive deep-inelastic scat- 
tering at HERA in the centre-of-mass (CM) frame of the photon dissociation (7* IP) system. The 
data support the conclusion reached in the analysis of the diffractive structure function F?^ 1 ' 
that, at low Q 2 , the momentum of the diffractive exchange is carried largely by hard gluons. 
Thus, significant transverse momentum and energy flow are produced by the hard subprocess 
in the boson-gluon fusion mechanism, and also arise from the enhanced soft gluon radiation 
associated with the gluon from the diffractive exchange. In the photon dissociation picture, the 
same data can correspondingly be interpreted as evidence for a significant contribution from 
events with photon fluctuation into Fock states with one or more gluons. 

The data have been compared with the 7*p data of the EMC collaboration, with the "f*p 
CM energy (W) of the EMC inclusive DIS data similar to the j*P CM energy (M x ) of the HI 
diffractive data. It is seen that additional transverse momentum is produced in diffractive scat- 
tering compared to the j*p case. A striking contrast is observed in the seagull plot, indicating an 
approximate symmetry between the target and current hemispheres in diffractive DIS, as would 
be expected in the hard-gluon picture of the pomeron and in the photon dissociation picture, 
whereas radiation is significantly suppressed in the proton-remnant region in non-diffr active DIS. 

The features of hadron production are well reproduced by a model featuring a factorisable 
pomeron flux in the proton if the partonic structure of the pomeron is dominated by hard gluons 
at the starting scale of Qq = 3 GeV 2 used in a DGLAP analysis of the diffractive structure 
function F^ 3 \ In contrast, a model with a pomeron consisting, at low Q 2 , primarily of quarks 
fails to describe the data in several respects. With respect to the 7*iP axis, the quark-based 
pomeron produces too little energy flow and particle production in the central region, but too 
much at large |r/| and \xp\. In the pj 2 and "seagull" distributions, this model produces too little 
transverse momentum. A model based on soft colour interactions (SCI) also gives an acceptable 
description of the data. 
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Figure 6: Seagull plot showing HI y*F> data in the ~/*]P CM frame in the kinematic region 7.5 < Q 2 < 100 GeV 2 , 
0.05 < y < 0.6, x P < 0.025, \t\ < lGeV 2 , My < 1.6 GeV and 8 < M x < 18 GeV, together with EMC fip DIS 
data in the 7*p CM frame. Positive xf corresponds to the direction of motion of the incoming photon. The 
points are plotted at the centre of the bin in the horizontal coordinate, and the horizontal bars indicate the width 
of the bin. Also shown are the predictions of the following Monte Carlo models for the HI data: RAPGAP with 
the hard-gluon pomeron structure taken from fit 3 in |] (RG-F° (fit 3)); RAPGAP with the pomeron structure 
containing only quarks at the starting scale, taken from fit 1 in Q (RG-Fjp (fit 1)); and the soft colour interaction 
model as implemented in LEPTO 6.5. 
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